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Let us consider the following operator (quantum Hamiltonian) in M. d with 

<N d = 3 

> 

(1.1) H= H AV = ((/)D-/\)-cr) 2 - \/(x) 

where /4, V are real- valued functions and V has a Coulomb- like singularity 
at or has several such singularities and is smooth and decays as Coulomb 

t-h or better at infinity 1 . 

Let A e Jff 1 . Then operator H is self-adjoint in «i? 2 (M 3 ,C 2 ). We are 
interested in Tr H A V (the sum of all negative eigenvalues of this operator). 

'X Let 

(1.2) E* = inf EM), 

(1.3) E(A) := (Tr H Ay + K^h" 2 J \dA\ 2 cbc) 
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1 In [15] we assumed that V is a smooth function. 
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with OA = (diAj) a matrix. 

This paper is the second step to the recovering sharper asymptotics 
of the ground state energy for atoms and molecules in the self-generated 
magnetic fields. 

Let Xj e M 3 (J = 1, ... , M, where M is fixed) be singularities ("nuclei"). 
We assume that 



(1.6) \D a W\ < C a z j{ £ j( x ) + !) \U X )) H v « : l«l < 2 



but at first stages we will use some weaker assumptions. Later we assume 
that V(x) decays at infinity sufficiently fast. 

In this paper we assume that k e (0, k*\ where < k* is a small constant. 
As k = we set A = and consider Tr~ H A V ; then our results will not be 
new. 

2 Estimates of the minimizer 

Let us consider a Hamiltonian with potential V and let A be a minimizing 
expression (1.3) magnetic field. We say that A is a minimizer and in the 
framework of our problems we will prove it existence. 

2.1 Preliminary analysis 

First we start from the roughest possible estimate: 

Proposition 2.1. Let k < k*. Then the near-minimizer A satisfies 



(1.4) 





X; 



zj > 0, zi + ... +z M x 1, 



l<j<M 



(2.1) 




and 



(2.2) 



\\dA\\ < Ck* 
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Proof. Definitely (2.1)-(2.2) follow from the results of [EFS3] but we give 
an independent easier proof based on [15] . 



(i) First, let us pick up A = and consider Tr(^E(0)^J with cut-offs 
^(x) = ip((x - xj)/£) where ip G ^q°{B(0, 1)) and equals 1 in 6(0, §). Here 
and below E (r) = 0(r — Ha ,y) is a spectral projector of H. 

Then 

(2.3) | Tr(^H AU (0)^) | < Ch' 2 as £ = L ■= h 2 . 

On the other hand, contribution of B(x, £) with £(x) = | miny |x — X/| > I* 
to the Weyl error does not exceed Cp 2 ^ 1 = Cp^lh^ 1 where h = h/p£ in the 
rescaling; so after summation over i > I* we get 0(h~ 2 ) provided p 2 < Ci^ 1 . 
Therefore we arrive to the following rather easy inequality: 

(2.4) | J (tr ei (x, x, 0) - Weyl^x)) dx\ < Ch~ 2 . 

This is what rescaling method gives us without careful study of singularity. 

(ii) On the other hand, consider A ^ 0. Let us prove first that 

(2.5) Tr(^W^) > Ch' 2 - Ch' 2 J \dA\ 2 dx as £ = 4- 



Rescaling x4x/f and r i— > t/£ and therefore h i— )■ /?£~5 x 1 and /4 i— > /4£5 
(because singularity is Coulomb- like), we arrive to the same problem with 
the same k (in contrast to section 4 of [15] where k h- >• k£ because of different 
scale in r and h) and with £ = h = 1. 

However this estimate follows from the proof in section 3 of [ES3] of 
Lemma 2.1, namely from (3.19)-(3.22) with Z = d = 1. 

(iii) Consider now ipe as in (i) with £ > £*. Then according to theorem 4.1 
of [15] rescaled 

(2.6) Tr (ipeH A]V ipi) > -Cp z £h- X - Ch~ 2 [ \dA\ 2 dx. 

Jb(x,2£/3) 

Really, rescaling of the first part is a standard one and in the second part we 
should have in the front of the integral a coefficient K~ 1 h~ 2 p 2 x p~ 2 £(h/ p£)~ 2 
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where factor p 2 comes from the scaling of the spectral parameter, factor p~ 2 
comes from the scaling of the magnitude of A, factor £ = £ 3 x £~ 2 comes from 
the scaling of dx and d respectively, and h/(p£) is a semiclassical parameter 
after rescaling. So, we acquire a factor p 2 £ < C. 

Then 

(2.7) j(tre 1 {x,x,0)-\Ney\ 1 {x))dx>-Ch- 2 -Ch- 2 j \dA\ 2 dx 

and adding magnetic field energy we find out that the left-hand expression of 
(2.1) is greater than the same expression with A = plus (C — k -1 )/7~ 2 ||<9/4|| 2 
minus Ch~ 2 which implies (2.1) and (2.2) as A is supposed to be a near- 
minimizer. □ 

Remark 2.2. We are a bit ambivalent about convergence of J Weyl-^x) dx 
at infinity, as for Coulomb potential it diverges. In this case however we can 
either replace Ha,v by Ha,v + V with a small parameter rj > or consider 
the left-hand expression of (2.1) plus magnetic field energy as an object to 
minimize. 



2.2 Rough estimate to a minimizer. I 

Let us repeat arguments of subsection 1.3 of [15]. Let us consider equation 
for an minimizer A as in (1.13) of [15]: 



(2.8) AA = -2Kh 2 J2(°j°k{hD k - A k ) x + u k <Jj{hD k - A k ) y )e{x,y, r)\ 



If we scale with the scale x i— > x/£, r i— > r/p 2 , h i— > H = h/(p£) then (2.£ 
would become 

(2.9) AA = 

- 2Kp 3 £h 2 ^2(ojC> k (hD k - p^A k ) x + a k CTj{hD k - p~ 1 A k ) y )e{x,y,r)\ y=x 

k 

and since so far p 2 £ = 1 we arrive to 

(2.10) AA = 

- 2nph 2 ^2(ajG k (hD k - p^A k ) x + u k Vj(hD k - p~ 1 A k ) y )e(x, y,r)\ y = x . 
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(i) Plugging for u = E(0)f and repeating arguments of [15] we conclude that 
in the rescaled coordinates 

(2.11) \\HD x u\\ < \\{(hD x - p~ l A x ) ■ cr)u\\ + Cp _1 ||>A|| 6 ||u||3 < 

\\{{hD x - p~ x A x ) ■ tr)u\\ + Cp~ x K~* ||/A|| 6 ||i/||3 ■ \\hD x u\\^ < 

\\{{hD x - p-M x ) • <t)u\\ + \\\hD x u\\ + Cip-^WAW'^WuW 

where ||^||6 calculated in the rescaled coordinates is ^™ 3 / 5 ||>4||6,orig (where 
subscript "orig" means that the norm is calculated in the original coordinates) 
which does not exceed C£~2 ||<9/4|| or j g < C£~2K2 due to (2.2) and therefore 

(2.12) \\hD x u\\ < C(l + p~h-h- l *K^\\f\\. 

Continuing arguments of section 1.3 of [15] we conclude that in the rescaled 
coordinates 

(2.13) hWAdAlUB^) + IIA^lUe^) < Kf\ 

k ■.= c«ft _1 (i + p- x h-*tr*K*y. 

Then either 

(2-14) II^IL. fl (x,f) + II^L, e( x,|) < c/cri 

or 

(2.15) + II^H^g^D < c\\dA\\ = c\\dA\\ orig r- 2 < c K k-l 

where in the rescaled coordinates 

(2.16) ||fi||* := sup |fi(x) - B(y)\ ■ |x - + | log |x - y\\y\ 

x,y 

In the latter case (2.15) we have in the original coordinates 
(2-17) l|M|L.fl(x.|£) < C ^~ l 

and we are rather happy because then the effective intensity of the magnetic 
field in B(x, £) is p~ 1 £\\dA\\ 00iB (x,e) < C^ 5 if we take p = £~2. 
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In the former case (2.14) let us consider (still in the rescaled coordinates) 
(3(x) = \dA(x)\£2. Then (3(x) has the same magnitude (3(y) in 7-vicinity of 
y with 7 = ef3(y)K^ 1 \ log(/3(y)K _1 )| _1 (or 71 = e, whatever is smaller). But 
then in the rescaled coordinates 

r 1 /3 2 (/3/<- 1 |log(/3/<- 1 )|- 1 ) 3 < C\\dA\\ 2 < C\\dA\\ 2 or - lg r x < CkW 1 
and then 

f3 < Ck&K*\ log(/3/<- 1 )|f 

which implies 

(2.18) (3 < Ch-l\\ogh\l, h=ht^ 
(as 7 x 1 the same arguments lead us to (2.17)). 

Therefore in the first round of our estimates we arrive to the estimates in 
the rescaled coordinates 

(2.19) \dA\<(3rK (5 := Ch-I- 5 

where we just estimated | \ogh\ by h~ Sl ; below we increase S if needed but it 
still remains an arbitrarily small exponent. 

(ii) In the second round we do not invoke \\A\\q but rather ||^||oo,B(y,7) < 
C/3£ _ ^7 where we consider a ball of radius 7 < 1 in the rescaled coordinates 
(and subtract a constant from A if needed), resulting in 

||A^||oo,e(x, 7 ) < CkH- 1 p(1 + /3-fi-^p- 1 ) 4 . 

Let us increase p to p' = (/3hl~i)^ = C£~^(/3h)^ = C£~^h~^~ 5 > £~i and 
use 7 = hjp'l = hr ~ 5 < 1. Then we arrive to 

(2.20) ||A/l|| 00 , s(x , 7) < KrK K := Ck/T^. 

Repeating arguments of the first rounds we conclude that either (2.17) holds 
or 

(2.21) \dA\<(3rK (J ■= kI- 5 < Ch-%- 5 ; 

then rescaled magnetic field is 0(/3£~^ / p) = 0(h~ 21 / 25 ~ 5 ). Here we returned 
to the natural scale (£,p) with p = £~2. 

(iii) One can also run third etc rounds, using partially arguments of sub- 
section 2.1 of [15]; then the rescaled magnetic field is 0{h~ s ). However to 
prove that the rescaled magnetic field 0(1) we need to modify them, and 
we do it in the next subsection. 
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2.3 Rough estimate. II 

In this step we repeat arguments of subsection 2.1 of [15] but we have a 
problem: we cannot use p = H^Hoo as we have domains X r = {x : £(x) > r} 
rather than the whole space. So we get the following analogue of (2.19) of 
[15] in the rescaled coordinates: 

(2.22) WAAW^^ + hWAdAW^K 

Ck P (ji + /r 1 /^- 1 ^-! ||&4|||* {fl{x>1)) ) 

with p = max(yU, 1), and p = p~ 1 |cM|oo I B(x,i)- But then 

(2.23) h^p-^dAW^^ < eh^ p - l \\dA\y e{B{xA)) + C K p+C K . 

Obviously in the right-hand expression we can replace p = P ~ 1 \dA\ OQi B( x ,i) 
by any other norm, in particular by 2 -norm 

P = P _:L ||94||b(x,1) = p~ l l~^\\dA\\ OT ^ 

which would be less than Cnh. 

Let u(r) = sup x:i , x ^ >r f(x) where f(x) is the left-hand expression of (2.12) 
calculated for given x in the rescaled coordinates. Then (2.23) implies that 

v{r) < -i/(-r) + C«s 

which in turn implies that 

v{r)<^v{2-"r) + 2C, n > 1, 

and therefore 

z/(r) < ACk* + 4 sup f(x) < Qn^ 

C h 2 <l(x)<2C h 2 

due to the rough estimate (because h~ x 1 as £(x) x h 2 ). Then going to the 
original coordinates we arrive to estimates below: 

Proposition 2.3. Let k < k* , p = c£~i. Let A be a minimizer. Then for 
£(x) >£ Uf = h 2 (2.17) holds and also 

(2.24) \d 2 A{x)-d 2 A{y)\ < OciH|x - y\H e l 2 £Z e ' 2 < 6 < 1, 
and 

(2.25) \dA(x) - dA(y)\ < CK*f*\x -y|(l + | log |x -y||). 
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Remark 2.4- (i) So far we used only assumption that 

(2.26) \d a V\ < Cp 2 r lal Va : \a\ < 2 
with p = but even this was excessive. 

(ii) In this framework however we cannot prove better estimates as (2.17) 
always remains a valid alternative even if p <C £~ 2 . 

(iii) Originally we need an assumption (2.4) of [15] | V\ > eo, but for d > 3 
one can easily get rid off it by rescaling technique; see also corollary 2.3(h). 

Consider now zone {x : £(x) < 

Proposition 2.5. Let k < k* , p < ci~^ . Let A be a minimizer. Then 
\dA\< Ch- 3 as £{x) < 4- 

Proof. Proof is standard, based on rescaling (then h = 1) and equation (2.8) 
for A. We leave details to the reader. □ 

Let us slightly improve estimate to A. We already know that |cM(x)| < 

3 

Co/3 with p = £~z and using a standard rescaling technique we conclude 
that 

(2.27) \AA\ < Ck P 2 /3 + Cup^t 1 

5 

which does not exceed Ck£~z which implies 
Proposition 2.6. In our framework 
(1) As £{x) > h 2 

(2.28) \A\<CkTK \dA\<CKT3, 

(2.29) \dA(x) - dA{y)\ < Qk£-\- q \x - y\ e as \x - y\ < -£{x) 

for any 9 G (0, 1) 

(ii) as £(x) < h 2 these estimates hold with £(x) replaced by h 2 . 
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Here in comparison with old estimates we replaced factor nh by k which 
is an advantage. 

Consider now zone {£ > max(a, 1)} and assume that 

(2.30) p < CC V as £ > ^ 

with v > 1. Then if also /3 = 0(£~ Ul ) as I > 1 the right hand expression 
of (2.27) does not exceed Ck^ -3 " -1 + £-^~ 2 ") and therefore we almost 
upgrade estimate to f3 to 0(£~ 3v + £~ Vl ~ 2u+1 ) and repeating these arguments 
sufficiently many times to 0(£~ 3v ). However, there are obstacles: first, as 
v > 1 we get 

Aj = J2° l jJ x -xn,]- 1 + oir 1 - 6 ) 

m 

with constant <Xj m ; however assumption V • A = implies 0Cj m = and we 
pass this obstacle. The second obstacle 

A J = a jl<,™( X k ~ *k,m)\x - XmT 3 + 0(r 2 ) 

k,m 

with constant <Xjk, m we cannot pass as assumption V • A = implies only 
that modulo gradient A = J2 m x with constant vectors |3 m and 

one cannot pass this obstacle. 

Therefore we upgrade (2.28)-(2.29) there: 

Proposition 2.7. In our framework assume additionally that (2.30) holds. 
Then as v > | 

(2.31) \A\ < CKt 2 , \dA\ < Ck£~ 3 , 

(2.32) \dA{x) - 8A{y)\ < C 6 Kr 3 - e \x - y\ e as \x - y\ < ht(x) 
as £{x) > 1 (for all e (0,1);. 

Remark 2.8. (i) In application we are interested in v = 2; 

(ii) We cannot improve (2.31)-(2.32) no matter how fast p decays. 

3 Tauberian theory 

Recall that the standard Tauberian theory results in the remainder estimate 
0(h~ 2 ). Really, as the rescaled magnetic field intensity is no more than 
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Ck2, contribution of B(x,£(x)) to the Tauberian error does not exceed 
Cp 2 x frr 1 = Cp 3 £h~ 1 which as p x translates into Ci~^h~ x and 
summation over {x : £(x) > £* = h 2 } results in Ch~ 2 . On the other hand, 
contribution of {x : £(x) < £* = h 2 } into asymptotics does not exceed 
C/r 3 ^ 1 = C/r 2 as h = 1. 

However now we can unleash arguments of [IS]. Recall that we are 
looking at 

(3.1) M^H AV ^) = Tr((p 1 H AV (f) 1 ) + Tr(0 2 H^0 2 ) 

where i\) 2 = <p\ + 2 , supp0i C {x, |x| < 2/?}, supp</> 2 C {x, R < \x\ < a} 
and we compare it with the same expression calculated for H AV o with 
V° = z\x\~ l . Here we assume that 

(3.2) a<l, zxl 
and 

(3.3) \D a (V - V°)\ < Coa-H-W Va : |a| < 3. 

The latter assumption is too restrictive and could be weaken. Then 

(3.4) Tr(MH AV -H Av0 ) ( p 2 ) = 

J (Weyli(x) - Weyl°(x)) 2 (x) dx + 0{R~h- 1 ) 

where WeylJ and Weyl are calculated for operator with potential V°. Really, 
we prove this for each operator H A y and H A V o separately 2 . 

On the other hand, considering \A = V°(l — () + V( = V° + W( and 
following [IS] we can rewrite the similar expression albeit for c/> 2 = 1 as 

(3.5) Tr f We(-H Ay c)d( 

Jo 

and applying the semiclassical approximation (under temporary assumption 
that W is supported in {x : |x| < 4/?}) one can prove that as <f>i = 1 

(3.6) Jr(4> 1 {H- AV -H- v0 )4>,) = 

J (Weyl^x) - Weyl?(x)) 2 (x) dx + O^ 1 /?/?- 2 ). 

2 Sure, such formula requires two-term expression but one can verify easily that the 
second term is 0. 
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Really, contribution of ball B(x, £(x)) does not exceed Ca~ l h~ 2 = Ca~ 1 £(x)h~ 2 
and summation with respect to partition as £(x) < R returns Ca _1 /?/)~ 2 ); 
meanwhile contribution of {x : £(x) < £*} does not exceed Ca~ l h~ 2 = Ca~ x 
as there h = 1. 

One can get rid off the temporary assumption and take 0i supported in 
{x : £{x) < 2R} instead. 
Therefore we arrive to 

Proposition 3.1. Under assumption (3.3) 

(3.7) Tr(VW, v - H- v0 )^) = 

J (Weyl^x) - Weyl?(x)) ^ 2 (x) dx + 0(a~h-l) 

Really, a~^h~t is R^^h^ 1 + a^Rh^ 2 optimized by R x /?* := (a/?)i; as 
h 2 < a we note that /? 2 < /?* < a. 

Corollary 3.2. (%) vis M = 1 equality (3.7) remains valid with -0 = 1 and 
a = 1. 

(wj y4s M > 2 an<i a > h 2 equality (3.7) becomes 

(3.8) Tr(iK>£ v - H- v0 )^) = 

/ (Weyl^x) - Weyl?(x)) ^ 2 (x) dx + 0((a"3 + l)/7~5) 

where we reset case a > 1 to a = 1. 

Remark 3.3. One can apply much more advanced arguments of [13] or 
section 12.5 of [14]. Unfortunately using these arguments so far I was not 
able to improve the above results unless k < 1. More precisely, I proved 
estimate 0(nh~ s ~i + h^ 1 ) as a = 1 (or even o(h~ l ) as a ^> 1, k = o(ti3 +s ) 
and some assumptions of global nature are fulfilled). However as I still hope 
to improve these results, I am not including them here. 
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4 Single singularity 
4.1 Coulomb potential 

Consider now exactly Coulomb potential: V = z|x| _1 . Then according to 
Theorem 2.4 of [EFS3] as h = 1, z = 1 and < k < k* 



(4.1) Jjmfinf (Tr(> R H^ R ) + i ^ |«94| 2 c/x^ 

~y Weyli(x)^(x) cfec^ =: 2z 2 S{zk). 
which according to Lemma 2.5 of [EFS3] coincides with 

(4.2) Mm 6nf ("Tr - (H A + 77) + ^ ^ |<94| 2 dx^ 



J \Ney\ 1 {H A + r],x)dx^j =2z 2 S{zk). 

Here cj) e ^{B(0, 1)), = 1 in 6(0, §), </> R = 0(x//?). Also due to scaling 
for z > one has a Scott coefficient 2z 2 S(nz). 

Proposition 4.1. ^4s < k < «' 

(4.3) S(k') < S(k) < S(k') + C/^kT 1 - k'- 1 ). 

Proof. Monotonicity of S(/t) is obvious. 

Let < k < k' < k" < k*. Then for any e > if R = R £ is large enough 
then the left-hand expression in (4.1) for k' (without inf and lim) is greater 
than S(k") — e + — /t //_1 )||9/4|| 2 ; also, if A is an almost minimizer there, 
it is less than S(k') + e. 

Therefore {k'- 1 - K"~ l )\\dA\\ 2 < \S(k") - S{k')\ + 2e. But then 

S(«) - e < S{k') + e + {k- 1 - k'-^WOAW 2 < 

S(k') + e + C( K - 1 - ac'" 1 )^" 1 - k"- 1 )- 1 ^^") - S(«')| + 2e) 

and therefore 

(4.4) (/sT 1 - k'- 1 )- 1 |S(«) - S(k')| < K -1 - «" -1 ) -1 |S(k') - 5(01 
which for «" = implies (4.3). □ 
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Remark 4- 2. Using global equation (2.8) we conclude that as 

(4.5) \d a A\ < C/tr 1 " 1 " 1 t>\, \a\ < 1, 

(4.6) \d a A\ < CkT^ £<1, \a\ < 1, 

(4.7) \\dA\\ 2 < Ck 2 . 
Then 

(4.8) S'(k)<C, |S(k(1 + 77)) - 5{r])\ < Ckt). 

4.2 Main theorem 

In the "atomic" case M = 1 we arrive instantly to 
Theorem 4.3. If M = 1, k < k* then 

(4.9) E* = J Weyl^x) dx + 2z 2 S(z/t)/T 2 + 0{h^). 

Proof. If /4 satisfies minimizer properties then in virtue of corollary 3.2 

(4.10) Tr~ Ha,v - J Weyl^x) dx = Tr H AV o - J Weyl?(x) dx 

4 

mod 0(/T 5 ) 

and adding magnetic energy and plugging either minimizer for V or for V° 

we get 

r 1 r x 

> 



(4.11) inf (V H A V - J Weyl^x) dx + J \dA\ 2 dx} 

inf (Tr H AV o - y Weyl?(x) dx+-^ J \dA\ 2 dx} ± Or 



Sure as V (and surely V°) are not sufficiently fast decaying at infinity the 
left (and for sure the right hand) expression in (4.10) should be regularized 
as in section 4. However for potential decaying fast enough (faster than 
|x|~ 2 ~ 5 ) regularization is not needed. 

For V° we have an exact expression which concludes the proof. □ 
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5 Several singularities 



Consider now "molecular" case M > 1. Then we need more delicate 
arguments. 

5.1 Decoupling of singularities 

Consider partition of unity 1 = X)o</<m ^/ wnere i>j is supported in la- 
vicinity of xy as j = 1, ... , m and ipo = in | a- vicinities of xy ( "near- nuclei" 
and "between- nuclei" partition elements). 

Estimate from above 

Then 

(5.1) Tr H A v = M1>jH*v1>j) 

0<j<m 

and to estimate E* from the above we impose an extra condition to A: 

(5.2) A = as £(x) > ^a. 

5 

Then in this framework we estimate 

(5.3) | Tr-(V>oH A ^ ) - j Weyl 1 (x)^ 2 (x) dx\ < Ch^a^. 

Proof is trivial by using ^-admissible partition and applying results of the 
theory without any magnetic field. 

So, to estimate E* from above 3 we just need to estimate from above 
minimum with respect to A satisfying (5.2) of expression 

(5.4) Tr^jH^j) - [ Weyl^x^x) dx + -L I \dA\ 2 dx. 

3 Modulo error in (3.8). 
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Estimate from below 

In this case we use the same partition of unity {ipj}j=o,i m an d estimate 

(5.5) Tr Hfy v > Tr^jH AV ^j) 

0<J<m 

with 

(5.6) V = V + 2h 2 J2( 9 ^) 2 

j 

and we also use decomposition 

(5.7) J\dA\ 2 dx= Ju]\dA\ 2 dx 



0<J<m ' 



with 



(5.8) uij(x) = 1 as £j(x) < ^a, K/'( x ) > ^ — Cq as £j(x) < ^a 

j = 1, ■■■ , m, 



(5.9) uj > e q as £(x) > -a. 

5 

So far q > is a constant but later it will be a small parameter. Then since 

(5.10) Tr-(^oH A v"^o) ~ J Weyl 1 (x)^ 2 (x) dx + J ^o\ dA \ 2 dx > 

(again proven by partition) in virtue of [15] we are left with the estimates 
from below for 

(5.11) Tr^jH^j) - J Weyl 1 (x)^?(x) dx + -L | w 2|^|2 ^ 

Remark 5.1. (i) Note that the error in Wey^ when we replace V there by 
V does not exceed C/j _1 (1 + a~5) which is less than error in (3.8). Here we 
can also assume that A satisfies (5.2); we need just to replace q by e ^ in 
(5.8)-(5.9). 

(ii) We also can further go down by replacing Tr~ (ipj H a, vipj) by Tr(ipjH^ v ,ipj). 
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(iii) Therefore we basically have the same object for both estimates albeit 
with marginally different potentials (V in the estimate from above and V 
in the estimate from below) and with a weight uj satisfying (5.8)-(5.9); in 
both cases w = 1 as £(x) < but in the estimate from above cu(x) grows 
to C and in the estimate from below co(x) decays to q as £(x) > |a and in 
both cases condition (5.2) could be imposed or skipped. 

(iv) From now on we consider a single singularity at and we skip index j. 
However if there was a single singularity from the beginning, all arguments 
of this and forthcoming subsections would be unnecessary. 

Scaling 

(i) We are done as z x 1 but as z C 1 4 we need a bit more fixing. The 
problem is that V x zi~ x only as |x| < za; otherwise V < a" 1 (where we 
assume that a < 1). To deal with this we apply in the zone {x : za < |x| < a} 
the same procedure as before and its contribution to the error will be Ch~ 1 a~? 
as p = a~2 here. Actually we also need to keep |x| > z -1 /? 2 ; so we assume 
that z" 1 /? 2 < za i.e. z > 3-2/7. 

Now scaling x h-> x' = x/za, multiplying H a V by a (and therefore also 
multiplying A by 32, so A 1— > A' = a^A, h 1— > h' = ha~*z~ x \ then the magnetic 
energy becomes k _1 /7 _2 z J u)(x) 2 \d' A'\ 2 dx' where factors 3 _1 and az come 
from substitution A = a~?-A' and scaling respectively. We need to multiply 
it by 3 (as we multiplied an operator); plugging h~ 2 = h'~ 2 a~ 1 z~ 2 we get 
the same expression as before but with z' = 1, a' = 1 and h' = ha^^z^ 1 < 1 
and k' = kz instead of h and k. If we establish here an error 0(h'-3) the 
final error will be 0(a~ 1 h'~i) = 0(a~ih~izi). 

(ii) On the other hand, let z < a~^h. Recall, we assume that a > Coh 2 . 
Then we can apply the same arguments as before but with z = a~ih 
and we arrive to the same situation as before albeit with h' = 1, a' = 1, 
k' = na~2h and with z' = z/z. Then we have the trivial error estimate 

0(3-*) = 0(3-1/7-1). 

4 As z denotes zj we assume only that z\ + ... + zm x 1. 
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5.2 Main results 

Combining results of the previous subsections with proposition 2.7 we arrive 
to 

Theorem 5.2. If M > 2, k < k* and (2.30) holds with v>\ then 
(5.12) E* = J Weyl^x) dx + 2 ^ z 2 S{ Zj K)h- 2 + 0(/?i + /? 2 ) 



with 

(5.13) /?i 



/?~3 a >l 
a -3/)-5 h 2 <a<l 



and 

(5.14) /? 2 = k/T 



a 3 3 > \ \0gh\3, 

I log /? 2 /a| 1 /? 2 < a< |log/i|s ' 

Proof. To prove theorem we need to prove an estimate 

(5-15) ^II^IIW)<26} < CR 2 

where < b < a is a "cut-off" . On the other hand we know that 

1 

(5.16) _ M > = -^=0(1) 

and we need to recover the last factor in the definition of /?2. 

As a > 1 we can have a -3 because in virtue of (2.31) the square of the 
partial norm in (5.16) does not exceed Ca _3 /t 2 . 

On the other hand, as h 2 < /?* < a we can select b : /?* < b < a such 
that he partial norm in (5.16) does not exceed C| log(a/ /7 2 )| 1 • ||<9/4|| 2 . □ 

Remark 5.3. (i) As a < | log/?| we do not need assumption (2.30); 

(ii) In particular, as a > 1 and k < a 3 /?i remainder estimate is 0(h~^). 
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5.3 Problems and remarks 

Problem 5.4. (i) As k G [0, n*] with small enough n* does S(k) really 
depend on k or S(n) = 5(0) ? 

(ii) If S(k) really depends on k, what is asymptotic behavior of S(k) — S(0) 
as k — y +0: can one improve S(n) — S(0) = O(n) ? 

Any estimate better than O(n) would improve (with respect to k) remainder 
estimates in theorems 4.3 and 5.2. 

Problem 5.5. Improve (as a > 1 ) estimates in theorem 4-3 and 5.2 to 

those achieved in section 12.5 of [14] for k = (i.e. without self- generated 
magnetic field) . Namely there we were able to achieve 0(h~ 1 ) or even better, 
up to 0(/r 1+<5 ) 5 . 

(i) The best outcome would be the same estimate 0(/? _1 ) (or better 5 ) for 
all k G [0, k*]. 

(ii) Alternatively, we would like to see estimate 0(h~ 1 + K fJ, h~^); in particular 
we would get estimate 0(1) for k = 0(h 2 ^ 3f1 ^) with exponent [i as large as 
possible. 
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